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Abstract 

The Rashba interaction induced when inversion symmetry is broken in solids is a key interaction 
connecting spin and charge for realizing novel magnetoelectric cross-correlation effects. Here, we 
theoretically explore the optical properties of a bulk Rashba conductor by calculating the transport 
coefficients at finite frequencies. It is demonstrated that the combination of direct and inverse 
Edelstein effects leads to a softening of the plasma frequency for the electric field perpendicular to 
the Rashba field, resulting in a hyperbolic electromagnetic metamaterial. In the presence of magne¬ 
tization, a significant enhancement of anisotropic propagation (directional dichroism) is predicted 
because of interband transition edge singularity. Based on an effective Hamiltonian analysis, the 
dichroism is demonstrated to be driven by toroidal and quadratic moments of the magnetic Rashba 
system. The effective theory of the cross-correlation effects has the same mathematical structure 
as that of insulating multiferroics. 
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I. INTRODUCTION 


The spin—orbit interaction is a relativistic interaction connecting the spin and orbital 
motion of charged elementary particles. This interaction plays essential roles in conversion 
between electric signals and magnetic ones [1, 2], Of particular interest is the Rashba 
spin—orbit interaction—the one that arises when the inversion symmetry of the system is 
broken [3]. It is represented by the quantum mechanical Hamiltonian 




(i) 


where m is the electron mass, h is the Planck constant divided by 2n, a is a vector rep¬ 
resenting the direction and strength of the Rashba field and p and cr are the momentum 
operators for the electron and vector of Pauli matrices, respectively. The interaction leads to 
various cross-correlation effects in electromagnetism, i.e. a mixing of electric and magnetic 
fields. A straightforward consequence of the Rashba interaction in the dc limit is the spin 
polarization of a conduction electron induced when an electric current is applied (called the 
Rashba—Edelstein effect) [4], and the electric current generated when a magnetic field is 
applied (the inverse Edelstein effect) [5]. The Rashba-induced spin polarization has been 
predicted to be highly useful for magnetization switching by electric current [6, 7]. Magne¬ 
toelectric effects have been intensively studied in insulating materials (multiferroics), where 
magnetization and electric polarization are mutually coupled [8]. A conducting Rashba sys¬ 
tem is unique in this sense, since magnetoelectric effects couple magnetization and electric 
current instead of electric polarization; thus, the effect is expected to be smoothly integrated 
into conventional electronic devices. 

Although the Rashba interaction realized in two-dimensional electron gases in semicon¬ 
ductors has generally been weak, a strong Rashba interaction of a = 3 eVA was recently 
observed at surfaces of noble metals doped with Bi [9]. Metallic surfaces and interfaces have 
now become highly important in spintronics for magnetization control and spin-to-charge 
conversion with high efficiency [10, 11]. It has been observed that topological insulators, 
the extreme cases of Rashba systems where the quadratic term in the dispersion disappears, 
have intriguing features such as a coupled collective mode of spin and charge (spin plasmon) 
on the surface [12], Bulk BiTel, a semiconductor where carrier electrons are introduced by 
Bi vacancies [13], was recently observed to have a large Rashba splitting of a = 3.85 eVA[14]; 
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a large magneto-optical (Kerr) effect [15] and a giant spin-polarized photocurrent [16] were 
also discovered. 

The coupling of electric and magnetic degrees of freedom by the Rashba spin—orbit inter¬ 
action leads to intriguing optical responses, including anisotropic propagation (directional 
dichroism) and electromagnetic metamaterial behavior; metamaterials are materials that 
exhibit extraordinary electromagnetic characteristics such as negative refraction [17]. Neg¬ 
ative refraction was originally discussed by Veselago and was considered to occur if the 
electric permittivity and magnetic permeability were simultaneously negative [18]. It has 
in fact been achieved by designing resonators for electric and magnetic responses with the 
same resonant frequencies in GHz [19] and THz [20] ranges. For shorter wavelengths, small 
structures comparable to the wavelength are necessary, and thus, artificial structures would 
not be suitable. Instead, natural materials become promising, since their structures can be 
designed down to the atomic size. In fact, oxide superconductors with layered structures, 
such as La 2 - a; Sr a ,Cu 04 and B^S^CaCi^Os+s, behave as metallic materials in the direction 
parallel to the layers and as dielectric materials in the direction perpendicular to the layers 
in the infrared frequency range [21], Such materials with strong anisotropic electric prop¬ 
erties (called hyperbolic materials) exhibit negative refraction when an interface parallel 
to the metallic plane is subjected to an electromagnetic wave. These materials are also 
useful for filters and for enabling hyperlenses to go beyond the diffraction limit [21, 22], 
Tetradymites, such as B^Tea and B^Ses, have been observed to be a new class of natural 
materials exhibiting a a hyperbolic nature for visible light [23]. 

In this paper, we present a quantitative and consistent theoretical analysis of the op¬ 
tical properties of a bulk Rashba conductor by evaluating the correlation functions from 
a microscopic standpoint. First, we shall demonstrate that bulk Rashba systems such as 
BiTel indeed exhibit naturally significant hyperbolic metamaterial properties in the infrared 
regime because of spin—charge mixing. These unique behaviors of Rashba systems are not 
observed in topological insulators, for which the bare plasma frequency is not defined by the 
quadratic dispersion. Second, we demonstrate that in the presence of magnetization or a 
magnetic field, the system exhibits directional dichroism because of breaking of both spatial 
inversion and time reversal symmetries. The effect is strongly enhanced by a singularity at 
the interband transition edges, the stronger of which is predicted to occur in the infrared 
and visible light regime for BiTel. Our analysis establishes a bridge from the theoretical side 
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between the optical properties discussed mostly on phenomenological grounds and solid- 
state transport phenomena based on the microscopic quantum mechanical viewpoint. The 
qualitative equivalence of the dc transport and optical responses is noted experimentally as 
well as from the symmetry viewpoints [24, 25]. The equivalence is natural from our analysis, 
since both phenomena are governed by a same conductivity tensor evaluated at different 
frequencies. 


II. CIRCULAR DICHROISM DUE TO SPIN-CHARGE CONVERSION 

When an electric held is applied, the Rashba interaction generates electron spin polar¬ 
ization (the so-called Edelstein effect) [4] (Fig. 1 a) whose density reads 

se = keol x E , (2) 

where ct = — is a unit vector along the Rashba held and Ke is a coefficient. The induced 
magnetization is Me = —H^se, where 7 is the gyromagnetic ratio. The interaction also 
induces the inverse effect, where an electric current density is induced by an applied magnetic 
held B. namely 


j ie — ftiE<d x B, 


( 3 ) 


where Kje is a coefficient. 

These cross-correlation effects also affect the optical properties. The inverse Edelstein 
current (3) is written by using Faraday’s law as (in the momentum and angular frequency 
representation) j'ie = 77 ( 0 : x (q x E)]. The Edelstein effect (2) also leads to a magnetization 
current density, j’e = V x Me- The coefficients of the Edelstein effect and its inverse effect 
are related as 


kie = ih'yujKE, (4) 

since both the spin density (2) and current density (3) are written in terms of the uniform 
component of the spin—current correlation function \sj (th e superscripts ij denote the 
direction). We thus see that the sum of the two generated currents is jm + jE = ih^KE^q^dt • 
E) — a(q-E)] and that the conductivity tensor has off-diagonal elements linear in q. The off- 
diagonal elements result in a rotation of the incident electric held around a rotation axis (3 = 
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TrfKviq x d) by an angle proportional to \(3\. The dispersion reads q 2 = ^±?'/io/3cu (/i 0 is the 
magnetic permeability) for a circular polarization ±, and hence, the light velocity depends 
on both the incident direction (sign of q x d) and the circular polarization (directional 
circular dichroism). This result is natural from a symmetry viewpoint, since breaking of 
inversion symmetry alone leads to the so-called natural optical activity, which affects only 
circularly polarized light [26]. 


III. HYPERBOLIC METAMATERIAL RESULTING FROM DIRECT AND IN¬ 
VERSE EDELSTEIN EFFECTS 


Peculiar behaviors of unpolarized light arise if the diagonal components of the conductiv¬ 
ity tensor are modified by the Rashba interaction that occurs when the direct and inverse 
Edelstein effects are combined. In fact, the magnetization Mp induced by the Edelstein ef¬ 
fect generates an electric current due to the inverse Edelstein effect: j’ie-e = ^m(d x M E ) = 
—h'yn E«iE[d x (d x E )]. The resultant polarization P = Jie-e/( — ■ ioj) is opposite to the 
applied electric held when E T d; in other words, the mixing of spin and charge excitations 
reduces the plasma frequency (Fig. 1 a). Microscopically, this softening effect is represented 
in the conductivity tensor in the long-wavelength (q —>■ 0) limit: 

aij{q = 0, u) = — . — (Sijil + C(u)) -OiidijCiu )V (5) 

ui + irj m \ J 

where 


CV) = 


4ci 2 


d 3 k 


n P 


-ep 


7 k^k 


( 6 ) 


(27r) 3 (Hu + irj) 2 — 4 

is a function representing the magnitudes of the Rashba—Edelstein and inverse Edelstein 
effects (plotted in Fig. 1 b). Here, Sk = XL=± a fka is the electron spin polarization 
(cr = ± denotes spin), = |'YfcI, 'Yfc = k x ct, f ka = [e £feCT d fc s T ) + l]- 1 [ s the Fermi 
distribution function, (ks and T are the Boltzmann constant and temperature, respectively), 


e ka = — cryfc — ep (cf and kp are the Fermi energy and Fermi wave vector, respectively) 

and rj is a damping parameter. 

In terms of C(u), the Rashba—Edelstein coefficient is k e = — lj m denotes 

the imaginary part). The function IniC^o;) is Enite for a broad frequency range lower than 
the higher interband transition edge cu + , where oj± = 4o;f5(\/1 + 5 2 ± 5) are the transition 
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edges, with 5 = being a dimensionless Rashba strength parameter and cuf = ^ (see 
Supplementary Information). For BiTel, ojp, the plasma frequency oo p , and are in the 
infrared regime, whereas oo + is in the visible frequency. 

Choosing a. along the z-axis, the dielectric tensor £ij(q — 0, uo) = S t] + = 0, a;) 

reads 


£7.7. 1 


£xx £yy 1 


00 , 


oj(oj + irj) 
,2 


P = 


oo„ 


00(00 + irj) 


— (1 + C(oo)) = £ x , 


(7) 


where oo p = e 2 n e /£o‘m is the bare plasma frequency (n e is the electron density). As 
observed in equation (7), the plasma frequency in the xy plane is indeed reduced by the 
Rashba interaction to wr = oj p ^1 + RuC^cur). 

The frequency range between the bare plasma frequency 0 o p and the reduced plasma 
frequency cur is of particular interest, since the system is insulating in the direction of the 
Rashba field but metallic in the perpendicular direction (Fig. 2 a,b). Materials with such 
anisotropy are called hyperbolic because of their hyperbolic dispersion [21]. In fact, choosing 
the wavevector q in the xz plane, the dispersion relation is c 2 g 2 — oo 2 £ x = 0 for E y (E l denotes 
the i-component of the electric field) and 

n 2 n 2 , ,2 

^7 + Qx _ 

£ x £ z C 


for E x and E z . When both Re e x and Re are positive, the dispersion is elliptic; however, 
it is hyperbolic if one of these values is negative (Fig. 2 c,d). Note that in equation (8) the 
component q z is governed by £ x with positive real part, whereas q x is by with negative 
real part; thus, the wave propagation in the z(x)-direction is insulating (metallic). 

Hyperbolic materials exhibit the following peculiar optical properties. When the inter¬ 
face is parallel to the metallic plane, the incident light tends to be refracted in the negative 
direction, since the metallic nature reflects the light (Fig. 2 a). This behavior is observed 
in the dispersion curve (Fig. 2 c), which reveals that the ^-component of the group velocity 
is negative. The Poynting vector S also has a negative ^-component. Notably, negatively 
refracted light is focused along the insulating direction, as observed in the plot of the trans¬ 
mittance defined by the ratio of the flow of S as T = (S in is the incident Poynting 
vector) (Fig. 2 e,g). In fact, the angle of refraction, determined by tan 9$ = Jr?, does not 
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exceed a threshold value d“ ax = tan -1 \ c ) anc ^ the absolute value of the thresh¬ 

old angle is small, since no solution exists for small |Re q z \ (Fig. 2 c). For uj/u p = 0.78, 
0 ma, x _ _]^ 3 ° (Fig. 2 g), the incident light even for a large angle of incidence is concentrated 
in a direction within 13° from the direction perpendicular to the interface. This focusing 
effect is stronger for smaller frequencies close to cur, whereas the distribution is broader for 
larger frequencies (Fig. 2 c,g,i). 

In contrast, when the interface is perpendicular to the metallic plane (Fig. 2 b), the light 
is mostly reflected when the angle of incidence is small (Fig. 2 f), since no solution exists 
for small \q z \ (Fig. 2 d). For transmitted light, in contrast to the parallel case, the angle 
of refraction cannot be small, and thus, the refracted waves propagate close to the interface 
(Fig. 2 h,j). An interesting feature of the refracted wave is the backward wave, where the 
wave vector q and the Poynting vector S are in opposite directions when the propagation 
is perpendicular to the interface (Fig. 2 b,d) [22]. 

For BiTel, the maximum plasma frequency known so far is = 2.5 x 10 14 Hz (corre¬ 

sponding to a wavelength of 7.5/iin) for n e — 8 x 10 25 m -3 and ep = 0.2 eV [15]. Using 
a = 3.85 eVAand kp ~ 0.1 A -1 , we have 5 ~ 1 , and thus, C(u) ~ —0.4 for frequencies 
uj < (Up (Fig. 1 b), resulting in = 0.77 (cur = 1.9 x 10 14 Hz, corresponding to the 

wavelength of 9.8/mi). Thus, hyperbolic behavior arises in the infrared regime. 

For experimental realization, the penetration depth A of the system, defined by A = Im j g j , 
is also essential. The imaginary part of the relative permittivity is on the order of 0.1 
in the hyperbolic regime, and the imaginary part of q z calculated from equation (8) is 
Im[g z ] ~ -Im[e]. For (jj/oj v ~ 1, we have A ~ 12/mi, and the absorption of the light is 
neglected if the sample thickness is smaller. 


IV. DIRECTIONAL DICHROISM OF THE MAGNETIC RASHBA SYSTEM 

So far, we have considered a pure Rashba conductor where time reversal symmetry is 
maintained. We now discuss the case without time reversal symmetry, as realized in a thin 
him of a Rashba conductor attached to a ferromagnet with uniform magnetization or in the 
presence of an external magnetic held. We thus consider the sd exchange interaction 

H sd = JsdM ■ a , (9) 


7 



where M is a dimensionless magnetization vector and J sc i represents the strength of the 
exchange interaction . 

The conductivity tensor including the effect of M to the linear order is (the coefficients 
are defined in the Supplementary Information) 


M _ 
a ij 


e 2 n P J. 


e^sd 


m 

i 


-4 D(u){a-M) 


^ijk^-k 


2 k^a 


K OLi£jkl “ 1 “ &j€ikl)^kQl 


6pA/pC<J \ 

E 3 (u) 


(e 1 


(oj) [(at x M) ■ q] [Sij — Q^Gj] “h 


E'2 (cd) 


(d x M)iq _lj f (ax M)jQ x,i 


_M ±(a x q)j + (a x q)iMf 


( 10 ) 


where q± = q — d(d ■ q) and M 1 - = M — (.M ■ d)d. The hrst term of equation (10) 
is proportional to (d ■ M) and corresponds to the anomalous Hall effect induced by the 
magnetization parallel to the Rashba held, which leads to an enormous magneto-optical 
response (Kerr effect) [15]. Other terms having diagonal components linear in the wave 
vector q induce directional dichroism when d x M is finite. We consider the case where the 
interface is in the xy plane with the magnetization along the y -axis and the light incident to 
the xz plane. For the electric held in the y-direction, the dispersion relation including the 
effect of magnetization (equation ( 10 )) leads to 


u: lWp J sd M 

q - -yfe + X—#i 3 (q ■ A), 

c 2 c 2 k F s F 


( 11 ) 


when —— —— < 1, where E\ 3 = E\ + E 3 , q = q/q and A = (d x M) (M = M/\M\) 

ck F 6 f 

is a moment that governs the dichroism. In the broad frequency range below ui + , 1111^3 is 
on the order of unity, whereas it has singularities at the transition edges (cu = w-t) (Fig. 1 
b). More precisely, the dominant contribution to the function ImFh 3 (a;) behaves as Im(^§, 
and, since IrnC oc 9(u± — c o)y/oj± — u> (9(x) is a step function) close to the edge at u±, 

3 

IrnKjs diverges at the transition edges if rj — 0 as Im£^i 3 oc (uj± —u)~ 2 . Considering a finite 
damping, the directional dichroism is strongly enhanced by a factor of 


IrnK 


13 


V 

hu y 


( 12 ) 


at the transition edges. Let us dehne the non-reciprocity ratio as the ratio of the transmission 

For the system with thickness t, 
0.25 x 10 -3 in the 


rate for the angle 6 q and n — 6 q as t/d = Teq — &q 


T e q +E-e q ' 


/q uiivi /i g 

= tanh , where t~ x = Irru/. For 6 q = 45° and t = 10/iin, ?/ D 



infrared regime with Im£i 3 ~ 1 for J^M/e f = 0 .1, which is already large enough for 
experimental detection. At the edge u±, //d is strongly enhanced to be approximately 0.24 
if ///(heap) = 0.01. 

Mathematically, the edge singularity results from the dichroism arising from the cor¬ 
relation function linear in q, which is obtained by expansion with respect to the angular 
frequency. An additional condition for imposing finite A was necessary for the singularity in 
the Rashba case.The singularity does not exist in the coefficient D, and only a cusp structure 
exists at ta_ (Fig. 3 b) (the ‘singularity’ in Ref. [27] is a cusp.). In the presence of of many 
interband transitions, the singularity is still expected to be present, since each transition 
contributes independently to the singularity. 

We have demonstrated that the non-reciprocal optical property is determined by the 
mutual orientation of q and a vector A. This vector is known as the Rashba-induced 
effective vector potential that drives electron spin flow [28-30], and in this context, the 
directional dichroism in the magnetic Rashba system is a result of the Doppler shift of 
the light due to the flow induced by the vector potential [31]. The factor q ■ A essential 
for directional dichroism is consistent with the observation of dc transport [25], where the 
electric resistance in the presence of applied electric and magnetic fields (E and B) was 
observed to be anisotropic depending on q ■ (E x B) (q is the wave vector of electron), as 
observed by the replacement a. —> E and M —» B. 

The effective vector potential A is also interpreted as a toroidal moment (a vector that 
couples with the rotation of magnetic held) and has been discussed intensively in multifer- 
roics [32], In fact, the isotropic contributions of our dichroism result (10) is described by an 
effective Hamiltonian 


H d = -gA ■ (E x B) - JiQijEiBj , (13) 

where g and h are coefficients and Qij = M t 6ij + Mj&i (see Supplementary Information). 
The interaction (13) is the interaction Hamiltonian proposed to describe cross correlation 
effects in insulator multiferroics, and the vector A and Q tJ in this context are toroidal 
and quadrupole moments, respectively [32], The form of interaction between the toroidal 
moment and the electromagnetic fields in equation (13) suggests that the edge singularity is 
a common feature of toroidal moment-induced directional dichroism. In fact, the energies of 
lights having positive and negative q ■ A get shifted according to the interaction; this energy 
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difference results in directional dichroism if the interband transition matrix element depends 
on the energy. As a result, the dichroism occurs proportional to the energy derivative of 
the interband matrix element, and the edge singularity arises whenever the matrix element 
vanishes at the edge ui e faster than linear, i.e. proportional to (c o e — uj with 0 < £ < 1. 

In the real-space representation, the isotropic part of the current responsible for the 
dichroism (contribution of equation (10)) is jd — V x M ( \ + jj_'-P^ , where M^,i = — g(A x 
E)i — hQijEj and Pd,i = x B)i + hQijBj are the effective magnetization and electric 
polarization induced in metals by the toroidal and quadrupole moments, respectively. These 
expressions are the metallic analogues of those in multiferroics, and our study indicates that 
the cross-correlation physics discussed so far for insulators is universal and can be applied 
to metals without changing their mathematical structure. 
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Figure 1 
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Figure 1. a: Schematic figure showing the cross-correlation effects in the plane perpendicular to the 
Rashba field a. The Edelstein (E) effect generates spin polarization Sg from the applied electric 
field E, and the inverse Edelstein (IE) effect generates electric current j’ie-e from the induced 
magnetization M# in the direction opposite to the applied field, b: The real (dotted lines) and 
imaginary (solid lines) parts of the function C(u) are plotted as functions of w/wp for a = 0.5,1,1.5. 
The interband transition edges ui± and imaginary part of £13 governing directional dichroism in 
the presence of magnetization are shown for the case of a = 1. c: Real and imaginary parts of 
e x (u;) and e z (cu) as a function of uj/lo p for a = 1.0 and 7j/(huj p ) = 0.01. The shaded region between 
wr and uip is the hyperbolic regime. 
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Figure 2. Optical properties in the hyperbolic regime for two geometries, the interface perpendicular 
to the Rashba held a (a,c,e,g,i) and the Rashba held within the interface plane (b,d,f,h,j). a,b: 
Schematic illustration of the geometries. The incident wave vector, the wave vector in the Rashba 
medium and the Poynting vector are denoted by qi n , q and S, respectively. c,d: Dispersion curves 
for a hxed angular frequency and the relation between the wave vector q and the group velocity v g . 
It is seen that the group velocity is refracted to the negative direction in the case c, while the wave 
vector q has a negative Re q x component (backward wave) in the case d. e,f: Transmittance T 
(solid lines) and angle of refraction for the Poynting vector 9$ (dotted lines) plotted as functions of 
the angle of incidence 9 m for several values of u)(= w/w p ) = 0.78, 0.88, 0.98, corresponding to close 
to the two edges and the center of hyperbolic regime. g,h: Transmittance T plotted as function of 
9s- Focusing effect and its dependence on the frequency are schematically shown in i,j. Focusing 
effect is stronger for larger wavelength (shown by red in i and j). 
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Figure 3 



Figure 3. a: Schematic illustration showing the configuration of the Rashba field a, magnetization 
M and wave vector q required for the directional dichroism to appear. The transmittance depends 
on the orientation of q with respect to A = (a x M). b: The imaginary parts of C{uj) and 
(representing the dichroism) and the real part of D{ 00 ) (representing the strength of the Kerr effect). 
The strength of the dichroism has a singularity at the transition edges uj±. whereas the Kerr effect 
has only a cusp. 
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Supplementary Information 


In this Supplementary Information, we present details on the calculations, which is rele¬ 
vant to the discussion of the theory in the main text. 


I. MODEL HAMILTONIAN OF THE MAGNETIC RASHBA SYSTEM 

In this section, we present the microscopic model Hamiltonian. Let us consider a bulk 
Rashba-electron conductor coupled with the magnetization, where the Rashba spin—orbit 
interaction breaks the space inversion symmetry and the magnetic coupling breaks the time 
reversal symmetry. The Hamiltonian of the conduction electron is given by 

H = c\. [e fe - e F + q • (cr X fc) - J s dM ■ a] c k , 

k 

= - e F + 7fc,M ■ cr] c fc , (S.l) 

k 

where c\ = (cj^, cj,,) is the electron creation operator with wave vector k , the up- 
arrow(down-arrow)f (f) represents the spin component along the ^-axis in the spin space, 
= K 2 k 2 / (2m) is the energy dispersion of free electrons, ep is the Fermi energy, cr is a 
vector of Pauli spin matrices, M is the magnetization vector and J s d represents the strength 
of exchange interaction and 'jkM = Ik — M with 7*, = k x a. The energy dispersion of 
the Hamiltonian in equation (S.l) is given by 

e k,M = e k + C7 fc,M, (S-2) 

where a = ±1 correspond to the spin splitting upper (a = 1 ) and lower (cr = — 1 ) bands, 
respectively, and 

7 k,M = sja 2 kl - 2 J sd {cx x M) ■ k± + J 2 d M • M. (S.3) 

where k L = k — [k ■ 6t)6t with 6t = a/|a| being the unit vector in the direction of Rashba 
held. When a x M = 0 the energy dispersion is symmetric with respect to A 7 -direct ion in 
fc-space. On the other hand, when a. _L M ^ 0, the energy dispersion becomes asymmetric 
in the fcj_-direction. This asymmetry due to the magnetic coupling implies that the light 
absorption for electronic interband transitions depends on whether the k± is parallel or 
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antiparallel to a x M, which will lead to the directional dichroism for electromagnetic wave 
propagation in the Rashba system. 

The one-particle Green’s function (resolvent) is given by 


Gk,M( Z ) 


(z ~ e fc + e F - 7 fc ,M ■ <r) 1 


1 

2 


E 

<T=± 1 


1 + QlkM • & 
Z ~ e k,M + e F 


(S.4) 


where ik.M — Ik.M / Ik.M- The retarded and advanced Green’s function are given by re¬ 
spectively, M (e) = Gk m (c + ?’ 0 + ) and G^ M (e) = Gk,M(c ~ * 0 + ), where 0 + is a positive 
infinitesimal. In this paper, we consider a clean limit and use the following relation in the 
latter calculations: 


Gk.M^G) ~ Gk_ M ( e ) ~ m ^2 d( e ~ e k,M + 6 f )(1 + (71k ,m ■ o), 

cr=d=l 


where S(x) is the delta function. 


(S.5) 


II. CHARGE AND SPIN RESPONSE 


Let us evaluate current and spin response induced by an external electromagnetic held 
based on the linear response theory. The perturbative Hamiltonian representing the linear 
coupling between current and spin of electrons and an external electromagnetic held is given 
by 


Hint(t) = - f dr | j(r) ■ A(r,t ) - frys(r) • B(r,t )}, 


(S.6) 


where j(r) is the current density operator, s(r) is the spin density operator, A(r,t) is 

d 

a vector potential yielding the electric held E(r,t) = —— A(r,t ) and the magnetic held 
B{r , t) = V x A(r,t), and y(= e/(2 m)) is gyromagnetic ratio. In the following calculations, 
we evaluate the expectation values of j(r) and s(r) in the Fourier space, whose components 
are given by 


jq = (v + a x o/h)c k+ , (S.7) 

k 

^ = crc fc+, (s.8) 

k 
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where v = hk/m is the conventional velocity, a x cr/h is the anomalous velocity due to the 
RSOI, which is related to the spin density operator and k± = k ± q/2. Fourier components 
of the electromagnetic held are also given by using A(q, co) [///. dt f dre~ ig ' r+iUJt A(r, t)^j 
as E(q,co) = icoA(q.co) and B(q,co) = iq x A(q,u), respectively. 

The expectation values of j q and s q are written on the linear response of E(q,oo) and 
B(q, co) as [SI] 


( jq,i)u 
(Sq,i)ui 


-ie 


co + i0+ L 


n P 


m 


X*jj{q,u;M) - -^Sij E j (q,u) + eh'yx t j a ( 1 L UJ > M ) B j( 1 li UJ )i ( S - 9 ) 


JS\ 


ie 


CO + l0+^ 


xl J j (q,u-,M)E j (q,u) - M )Bj(q,u), 


(S.10) 


where e > 0 is the elementary charge and n e is an equilibrium electron density and Xjji 
X l j S , xl] and x l J s are current—current, current—spin, spin—current and spin—spin correlation 
functions, respectively. Here and hereafter, repeated indices generally imply summation 
(i,j = 1,2,3). These are written by using the Green’s function as 


/ °° de 

—tr [v i G k+ .M(e + )v j G k _ M {.£-)\ 

/ °o ( J f 

—tr lv i G k+) M{.e+)v j G k _ tM {.e-)\ 

K 

){q,u\M) = Y, f 

k 

/ °° de 

—tr [. 'J i G k+M (e + )a j G k _ M {e -)] 


11 

. 

SJ 


2tt i 

de 

- 1 

2ni 


00 de 

x—tr [, rJiG k+ , M (e+)vjG k _ tM (e-)\ 


x 13 

Ass 


(5.11) 

(5.12) 

(5.13) 

(5.14) 


where v = v + ot x cr/h, k± = k ± q/2 and e± = edzco/2. The Green’s function G'fc jvr(e) 
here is a path-ordered one [S2]. The lesser component of G k ,m(c) is given by [S4, S3] 


— /( e ) ^fc,M( e )) > 


(S.15) 


where /(e) = (1 + e e / T ) 1 is the Fermi distribution function with T being the temperature. 

The key to understanding the optical magnetoelectric properties is to evaluate the optical 
conductivity, <Jij(q,co; M), which can be related to the correlation functions as 


<7ij(q,w) = 


—ie~ 


oo + i0+ 


A 7 7 w tJ 

. JJ m . 


'"e r 
- Oi 


+ 


ehx 


oo + i0+ 


Xjs £ lkj + £iUXsj 


Qk + 


ih 2 X 2 

CO + i0+ 


Xss ^ikm^jlnQkQh 

(S.16) 
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where £ij k is the totally anti-symmetric tensor with £123 = 1. In the main text, the optical 
conductivity aij(q,ui] M) is evaluated to the first order in q and M. Thus oj;M) 

may be expanded up to the linear in q and M. For ca; M ) and x^q.uj] M), we can 

put q = 0 and evaluate these up to the first order in M. Since the last term in equation 
(S.16) is the second order in q , this is negligible in the optical conductivity. 


III. CALCULATION OF OPTICAL CONDUCTIVITY TENSOR 


In this section, we present calculation results of correlation functions, Xjji Xj S i Xsj an( i 
Xij and the optical conductivity tensor, cr^, and give the explicit form of coefficients, C, D, 
Ei (i — 1,2, 3), of rrh in the main text. 

Let us first perform the e-integral in equations (S.11)-(S.14). Following the Langreth’s 
method [S4, S3], one can calculate the lesser component of G(e+)G(e_) as 


[G(e+)G(e_)] < = G<(e + )G A (e_) + G R (e + )G<(e_) 

= /(e + )(G A (e + ) - G R (e + ))G A (e_) + /(e_)G R (e + )(G A (e_) - G R (e_)). 

(S.17) 


Applying this rule to the correlation functions of equations (S.11)-(S.14) and using equation 
(S.5), we can carry out the e-integral due to the delta function. The next step is to expand 
correlation function with respect to q and M and subsequently to carry out the integration 
with respect to the rotation angle around the a-axis. Since these calculations are rather 
complicated, we here only show the results as 


ij n e ^ / N(-_|_ . Tl e Jsd]El w, -\/t\ 

Xjj — -<~Tca)<U - *-2 —D{u){at ■ M)e ijk a k 

JJ <YY 1 J nnr\ 


+ 


m 
Tie Jsd 


m ei 


m kpe p 


Ei(u)(a x M) ■ q8± + -E 2 (w) {(ax M)iqf + g A (« x M)j) 


h n, 


X'js — ~~ C ( UJ ) £ ijk & k + * 


. h n e J sd hu 


a m 


a m ei 


{Ni{u)a.j - a ■ N(u})S i: j ), 


ij ^ ,h n e JgdfiUJ . „ , . „ 1.J-/ \ r \ 

Xsj — --—C{u) £ ijk a k - 1- — — { a iNj(u) - OL ■ N[u)5ij ), 


a m 


a m ei 


», ± . h 2 n e J sd f\uj 

XL - . ( - 2 Sij) - i— — ^—£ ijk N k (uj) 


a m 


a* m ei 


(5.18) 

(5.19) 

(5.20) 

(5.21) 


6f 


where 8^ = 8^—a^aj with d = a/|a|, N(ui) = D(u)(M-a)a+a—E 3 (uj) (M — (M ■ a) a 


hu 
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and the frequency-dependent coefficients are defined as 



E 3 (oj) — 



( Sfc 87 k Sk \ 

\l k H k H k Hi )■ 


(5.22) 

(5.23) 

(5.24) 

(5.25) 

(5.26) 


with a = — being the dimensionless Rashba strength parameter, H k = (hu + i0 + ) 2 — 

tl /Cp 

and S k = Y,a=± a fka, n k = Y,a=± f k(T , s' k = J^a=± a fL > < = Ea=± fka with fka = 

df(e) 

f\ e k + <^7fc - e F ) and f ka = —— 

e=e k +a^ k -e ¥ 

Substituting equations (S.18), (S.19) and (S.20) into the equation (S.16), we obtain the 
optical conductivity as 
• o 

XC- Tl 

cTijiq,^; M) = — (hy(l + C(u)) - OLiOLjC (w)) + ih'yK E (uj)(q i a j - a^) + a™(q,u) 

(jJ I X U Til 

(S.27) 


, . . ehrig C(ui) . . 

where ^(cu) = 1 -is the Rashba—hdelstem coefficient and 

rna u 

Jx X 

a™(q,u) = - — ^-D(u)(d ■ M) £ ijk dk ~ + £ikidj)d k qi 


. e 2 n e 1 

\ - 

m k E e E u 


~-E 3 (u) [MT(a x q)j + (a x 


^i(w)(a x M) ■ qST + ~E 2 (u) {(ax M^qf + qf(a x M)j} 


(S.28) 


with q 1 = q — {q ■ a)a and = iVT — (iVf ■ a)a. The hrst term of equation (S.27) leads 
to the softening of plasma frequency in the directions perpendicular to ct. The second term 
represents the direct and inverse Edelstein effects, which leads to the natural optical activity 
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(see main text). The effect of the magnetization is contained in afj . The term proportional 
to (a ■ M) of equation (S.28) gives an anomalous Hall effect, which leads to a magneto¬ 
optical response (Kerr effect). Other terms including 6l x M or M 1 - and having diagonal 
components linear in the wave vector q of equation (S.28) induce directional dichroism. 


IV. WAVE PROPAGATION IN THE RASHBA SYTEM 


In this section, we consider the wave propagation in the Rashba system based on the 
optical conductivity of equation (S.27) and lead to the dispersion relation of equation (11) 
in the main text. 

When we choose a = (0, 0,1), M = (0, M, 0) and q = (q x , 0, q z ), the optical conductivity 
tensor can be written as 


a ij 


( + + 2 q x 0 ih^K E q x 

(1 + C)a D + i^^E 13 q x 0 


\ 


0 

-ihqn E q. x 


171 /CpCpCd 

0 u D 


(S.29) 


where E V2 — E\ + E 2 and E \ 3 = E\ + E 3 . Thus the dielectric tensor s l3 = Sij H- a tJ is 

£qUJ 

given by 


where 


/ 


£ ij 


£ x 0 


~0 


\ 


0 £y 0 

0 o £ z 


u‘J sd M , 

e x = £ x -£-- (Ei + E 2 )q x , 


Ey E X 


£* = 1 


CO 2 kpdp 
JsdM 
CO 2 kpdp 

uo{uo + iq) ’ 

,2 


(E\ + E 3 ) q x , 


£ =1 - — 

co(co + iq) 

0 = ^-0 
CO 


— (1 + C ), 


(S.30) 


(5.31) 

(5.32) 

(5.33) 

(5.34) 

(5.35) 
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with (5 = h^'K^p/ix■ In general, the electromagntic wave propagation in conducting media is 
described by the following equation: [S5] 

\c 2 (q 2 Sij - qiq-j ) - Ej = 0, (S.36) 


where c is the light velocity. Substituting equation (S.30) into the equation (S.36), we obtain 
the following characteristic equation yielding a dispersion relation for plane wave as 


r 2 n 2 , ,2 x 

c q z ~ u £ * 


0 


c\ql + ql)-co% 


~c 2 q x q z - Ho/3w 

0 


-c 2 q x q z 


ho An 


cW x 


= 0. 


(S.37) 


Thus for E y , which represents an ordinary wave, we have the dispersion relation 


c (q x + q z )-u £ y = 0. 


(S.38) 


Putting the wave vector on q = q(sin9 q , 0, cos 9 q ), we obtain 


uj luJ 2 J cd M . 

q - -vex + o — E ^q ■ A 

c 2 c l kpep 


(S.39) 


when —p— —— <C 1, where q = q/q and A = ol x M with M = M/\M\, which leads 

CKp 6p 

to equation (11) in the main text. Thus we see that the second term including the factor 
q ■ A of equation (S.39), which depends on the wave propagating direction q, leads to the 
directional dichroism. The difference in the absorption rate of positive (6 q ) and negative 
( 7 r — 9 q ) wave propagation is proportional to the imaginary part of if 13 (see main text). 


V. k- INTEGRALS OF COEFFICIENTS, C, D, E u E 2 AND E 3 


In this section, we show the results of ^-integrals of equations (S.22)-(S.26). 
Using 


5 -E 

k 


$k 

~ L‘ H k 


% 

H h 


(5.40) 

(5.41) 
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we can write equations (S.22)-(S.26) as 


Ciu) = C(0) - —e F (hu) 2 S 
n e 

4 Q? 

D(u) = - e 3 S 

n e 

O ^ 6p v ■v 


/ x LX CJh X -X / . ma . i 

EM = teTfc E + ir"‘ + iz c ^ 


2a 


3 9 dS 3 57V 

- 2S+ 2 iv -4^ + 4 w a7 


^ ci 1 e F / / ma 2 , \ 2a 9 / 1 3 5S 1 1 57V 1 

B2<w) = ~~ E V feS ‘ + ir"v + 7 fF [ s ~2 n + i w s: _ i w + s 


4n P hen 


-E 3 (o;) = +—e F hcn - TV + 

n P V 5cn 


where 


C(0) = - £P ^^. 

n e “ 7fc 


At zero temperature, we have 


where 


P 

A/T7 


n P = — 


y! y! <K e fc + CT7fc — £f) — 7^(1 + 2Aq), 


C(0) = 


k <j 

a 2 3n 2 


r»7T poo 


1 + 2A (5 A; F 47r 2 a Jo 


£p - / (75 / d/ch E crO (e fc + 0 " 7 fc - e F ) 


1 + 2A e 


A 5 = -a 2 i 1 + 


1 + a 2 


a 


tan 1 a 


At zero temperature, we write S in equation (S.40) and N in equation (S.41) as 


S = 


N = 


4vr 2 a J 0 

1 / 


47r 2 


/ OO 

dkk 

r poo 

d/0 sin 5 J dkk 2 ^ 


d(e fc + q~7fc ~ £f) 
(hu + 70+) 2 — 4q| 


S(e k + <77fc - e F ) 
(hu + 70+) 2 — 4q^ 


(5.42) 

(5.43) 

3 2 d 2 S 
8 U 

(s.44) 

2 &S\ 

du 2 J 

(5.45) 

(5.46) 

(5.47) 


(5.48) 

(5.49) 

(5.50) 

(5.51) 

(5.52) 
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Carrying out the integrals, we obtain the real and imaginary parts of S and N as 


ReS = 


1 1 
87 r 2 a 3 


U 


(A — u)L\ + 2 tan 1 


a 


A=±l 


Inns' = XS ' 


167 T a 3 oj 


A 


A=±l 


111 
167T a 3 u 


x < 


S + — S- 0 < LJ < o;_ 

— S— U- <lo < uJ + 

0 oJ + < u 


87 C cr A — oj 

A 


' 1 1 

— -h —— 0 < id < U— 


1 Ct 2 ^1 1 0( 2 

ImlV =-- > — =-- x < 

167tq; 3 ^ S\ 16n a 3 

a=±i A 


s + s _ 

1 

si 

0 


o ;_ < oj < oj. 


uj + < OJ 


(S.53) 


(S.54) 


(S.55) 


(S.56) 


where ujy = cu/(4cuf) with fouj = 6f is a dimensionless angular frequency normalized by the 
Fermi energy, 


L A = tan- 1 (^=) , 


Rx = 


■ tan 


-1 


m, 


VQ\ WO\J 

u 2 + 2Ab 2 a7 — b 2 

(A^P ’ 

S\ = \/a 2 ~ 2Ab 2 o; — ch 2 , 


Q A — 


(5.57) 

(5.58) 

(5.59) 

(5.60) 


and 


uj± = a (\/1 + b 2 ± b ) 


(S.61) 


are frequencies at the transition edges. The imaginary part of C(oj ) has cusps at cn± (see 
equation (S.65) bellow). Note that when Q\ < 0, the following replacement is understood 


1 (—Z=] ->■ 


VQIJ 


VQ; 


■ tan 


-1 


a 


VQ\ 


^ tanh " 1 fe) = 5 

1 , _i / d A 

—>■- , tanh 1 


y/-Qx\n 


V~Q a + <5 


V-Qa - <5 

V~Q\ + <5 


(S.62) 


V~Q> 


kV~Q> 


2V^Q~x 


In 


V-Qa - a | 

(S.63) 
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Substituting these results into equations (S.22) and (S.23), we obtain 


ReC=? a "‘ 


4 1 + 2A, 


U 


^ (A — u)L\ + 2a ; 2 tan 1 


a 


A=±l 


1 + 2 A„ ’ 


Im C = — - \ UjJ^XSx, 

8 1 + 2 A 5 ^ A ’ 


ReD = 


3 d 


-l 


16 1 + 2 A^ 


u 


yy (A — ou)L\ + 2 tan 1 


a 


A=±l 


Im D = — — -Y AS*, 
32 1 + 2Aa UJ 2 —^ 

A 

_ 2Ao,d _1 e F 3 d ' 2 

ReEi — -——— — h — - 


1 + 2Aa hjj 81 + 2A«5 


4a;tan 1 d + 2 y^(A — u)L\ H—d 2 


A — 2a; 


A — u 


(5.64) 

(5.65) 

(5.66) 

(5.67) 
A Ry 


— ^- O-A 


(A — 3a;)(1 + a 2 R\) 


— -a 


A 


a ; 2 + 2Ad 2 a; — d 2 


ImiA = 


Re Eh = 


37T d 


-2 


32 1 + 2Ac 


J2 XS - 


,-^ 2 Ev+v+ 2 


4 +^a 4 


\ ^ A 

or g or —an ^-3 
a + 


(5.68) 

(5.69) 


2 Ao,d 1 e F 3d 


-2 


1 + 2A 5 hu 81 + 2Ac 


2 a; tan 1 d + ^(A — o;)La + ^d 2 ^ ^—+Ai?A 


A — u 


1 _ 2 (A — 3a;) (1 + a 2 R\) 


— a 


A 


a ; 2 + 2Ad 2 a; — d 2 


Irn +9 = 


37T d 


-2 


32 1 + 2 A d 


ReE 3 = - 
Irri+3 = 


3 d 


-2 


u 1 + 2 Aq, 

37 t or 2 


E AS* - + 2 E # + + 2 < 2 “ 4 -a 2 - + E4 

.a A A A A . 

- 2 a; tan -1 d — (A — o;)La — d 2 —- —XR\ 

■A—^ ^ A — hi 


■2- 2 E 


A 


(5.70) 

(5.71) 

(5.72) 

(5.73) 


32 1 + 2A* V E a 

A 

Thus, from equations (S.69) and (S.73), we obtain the explicit form of IrriEi 3 = Im(Ei + E 3 ) 
as 

37T d -2 


Imi?i 3 = 


32 1 + 2Ac 


E - ^ 2 E ^ + j+ 2 “ 4 - a 2 - + E | 

A 


(S.74) 


We see that this function diverges towards the transition edges u>± because of the second and 
last terms. The dominant term at the edges is the last term. Thus the directional dichroism 
is strongly enhanced (see main text). Note that similar equations as (S.54) and (S.56) are 
understood in the equations (S.69), (S.71), (S.73) and (S.74). 
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VI. TOROIDAL AND QUADRUPOLE MOMENTS 


In this section, we briefly make an analogy between our microscopic results of af- in 
(S.28) and the conductivity induced by an effective toroidal and quadrupole moments. 

Let us consider an effective Hamiltonian proposed to describe cross correlation effects: 
|S6] 

H s = -gA ■ (E x B) — hQijEiBj (S.75) 

where g and h are coefficients, and A and Qij(= Qji) are a toroidal and a quadrupole 
moments, respectively. The induced current density is obtained from the derivative of the 
Hamiltonian with respect to the vector potential, A, as 

Jd,i = = crfjEj, (S.76) 

where the conductivity crfj is calculated as 

crfj = iguj [2 A ■ qSij - (,AiQj + )] + ihu ( Qu£ikj - SmQij) Qk- (S.77) 

In the real space representation, the induced current is written as 

HD 

J d = VxM d + —(STS) 
at 

where = —g(Ax E)i — hQijEj and P^i = g(Ax B)i+hQijBj are an effective magnetiza¬ 
tion and electric polarization induced by the toroidal and quadrupole moments, respectively. 
If we put respectively A and on A = a. x M and Q tJ = atiMj + we have 

crfj = 2 ig(a x M ) ■ - ig [(a x M) i q j + ^(a x M)*] 

+ ih [Mt(a x q)j + (a x q)iM.j] + ih [cq(M x q)i + (M x q)idj]. (S.79) 

Thus, except for coefficients, we can have derived the linear terms of afj with respect to q 
and M. Some differences between the two may come from the anisotropic property due to 
the Rashba field a. 
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